Exchange parameters of copper-based quasi-two-dimensional Heisenberg magnets 
measured using high magnetic fields and muon-spin rotation 
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Pulsed-field magnetization experiments (fields B of up to 85 T and temperatures T down to 0.4 K) 
are reported on nine organic Cu-based two-dimensional (2D) Heisenberg magnets. All compounds 
show a low-T magnetization that is concave as a function of B, with a sharp "elbow" transition to a 
constant value at a field B c . Monte-Carlo simulations including a finite interlayer exchange energy 
J± quantitatively reproduce the data; the concavity indicates the effective dimensionality and B c is 
an accurate measure of the in-plane exchange energy J. Using these values and Neel temperatures 
measured by muon-spin rotation, it is also possible to obtain a quantitative estimate of \J±/J\. In 
the light of these results, it is suggested that in magnets of the form [Cu(HF2)(pyz)2]X, where X 
is an anion, the sizes of J and J± are controlled by the tilting of the pyrazine (pyz) molecule with 
respect to the 2D planes. 

PACS numbers: 76.30-v, 75.10.Pq 
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I. INTRODUCTION 

Systems that can be described by the S = | 
dimensional (2D) square-lattice quantum Heisenberg an- 
tiferromagnet modet^ 2 -^ continue to attract considerable 
experimental 4 *^ and theoretical^^ attention. Recent 
impetus has been added to this field by suggestions that 
antifcrromagnctic fluctuations from S = h ions on a 
square lattice play a pivotal role in the mechanisms for 
superconductivity in the "high T c " cuprates^ ' 10 ' 11 ' 12 ' 13 ' 14 
and other correlated-electron systems 1 ^; moreover, 2D 
Heisenberg magnets have been suggested as possible test- 
beds for processes applicable to quantum computation 4 ^. 

Though long-range magnetic order cannot occur above 
T = in a true 2D Heisenberg system 2 -^, real materials 
that contain layers approximating to 2D Heisenberg sys- 
tems 2 -^ ' 17 ' 18 invariably possess interlayer coupling that 
can lead to a finite Neel temperature^ 7 -^. In this con- 
text, synthesis of organic complexes containing ions such 
as Cu 2+ , neutral bridging ligands^ and coordinating an- 
ion molecule s 17 ' 18 has proved fruitful in the production of 
a variety of one- and 2D magnetic system s 17 ' 20 ' 21 ' 22 . The 
current paper describes high-field magnetization mea- 
surements on nine Cu-based quasi-2D Heisenberg mag- 
nets that employ pyrazine (pyz) as a neutral bridging lig- 
anc ji7 i i8 i 23_ rp^g £ a ^ a s y 10w that the field (independent, 
low-tempcrature (T) magnetization M(B) shows a char- 
acteristic sharp "elbow" feature at the transition to sat- 
uration, with a concave curvature at lower B. Monte- 
Carlo evaluations of a 2D Heisenberg square lattice with 
an additional interlayer exchange coupling energy Jj_ re- 
produce the data quantitatively; the degree of concavity 
depends on the effective dimensionality of the system, 



whilst the field at which the "elbow" occurs is an accu- 
rate measure of the in-plane exchange energy J. Using 
these J values in conjunction with Neel temperatures de- 
duced from muon-spin rotation (/xSR), it is then possi- 
ble to gain a good estimate of the exchange anisotropy 
\Jj_/J\ for all of the magnets. 

Having established these findings using the whole 
range of compounds, we suggest that in magnets of the 
form [Cu(HF2)(pyz)2]X ) 18 ' 19 where X is an anion, J and 
J± may bc influenced by the tilting of the pyrazine (pyz) 
molecule with respect to the 2D planes. 



II. EXPERIMENTAL DETAILS 

The quasi-2D magnets studied in this work are listed 
in Tabic |IJ The samples arc produced in single or poly- 
crystalline form via aqueous chemical reaction between 
the appropriate CuX2 salts and stoichiometric amounts 
of the ligands; further details are given in Refsj 18 ' 23 ' 24 , 
where structural data derived from X-ray crystallogra- 
phy are also found. Some compounds had crystals large 
enough to permit measurements with a single, orientated 
sample (Table|T|. In other cases, the materials were either 
polycrystalline or the crystals too small for accurate ori- 
entation; therefore samples composed of many randomly- 
orientated microcrystals, effectively powders, were used. 
In addition to the characterization described in Rcfsi 18 ' 23 , 
sample (/-factors were measured 2 ^ using standard elec- 
tron paramagnetic resonance (EPR) 2 ^. 

The pulsed-field M experiments used a 1.5 mm bore, 
1.5 mm long, 1500-turn compensated-coil susccptometer, 
constructed from 50-gauge high-purity copper wir o 27 ' 28 . 
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Compound 


Be (T) 


a 


J(K) 


Tn (K) 


\J±/J\ 


[Cu(HF 2 )(pyz) 2 ]BF 4 


18.0 


2.13 ±0.01 


-6.3 


1.54 


9 x 10" 4 


[Cu(HF 2 )(pyz) 2 ]C104 B x 


19.1 


2.30 ±0.01 


-7.3 


1.94 


2 x 10" 3 


[Cu(HF 2 )(pyz) 2 ]C104 B ]{ 


20.9 


2.07 ±0.01 


-7.2 


1.94 


2 x 10~ 3 


[Cu(HF 2 )(pyz) 2 ]PF 6 


35.5 


2.11 ±0.01 


-12.4 


4.31 


1 x 10~ 2 


[Cu(HF 2 )(pyz) 2 ]SbF (i 


37.6 


2.11 ±0.01 


-13.3 


4.31 


9 x 10~ 3 


[Cu(HF 2 )(pyz) 2 ]AsF 6 


36.1 


2.16 ±0.02 


-12.9 


4.34 


1 x 10~ 2 


[Cu(pyo) 2 (pyz) 2 ](C10 4 ) 2 B ± 


20.8 


2.30 ±0.01 


-7.9 


< 2.0 


< 10" 3 


[Cu(pyo) 2 (pyz) 2 ](C104) 2 S|| 


22.2 


2.07 ±0.01 


-7.6 


< 2.0 


< 10 -3 


CuF 2 (pyz)(H 2 0) 2 B ± 


33.1 


2.09 ±0.02 


-11.6 


2.54 


3 x 10~ 4 


CuF 2 (pyz)(H 2 0) 2 B n 


28.8 


2.42 ±0.01 


-11.5 


2.54 


4 x 10~ 4 


Cu(pyz) 2 (Re0 4 ) 2 


42.7 


2.10 ±0.05 


-14.9 


4.2 


3 x 10~ 3 


Cu(pyz) 2 (H 2 0) 2 Cr 2 7 


13.3 


2.13 ±0.01 


-4.7 


< 1.6 


< 1 x 10 -2 



TABLE I: The quasi-2D magnets studied in this work, along with their saturation fields B c and g-factors (here pyz is pyrazine, 
pyo is pyridine-N-oxide) . Data for oriented single crystals are indicated by B\ | (B parallel to 2D layers) and B± (B perpendicular 
to 2D layers); other data are for powders. In the latter cases, an average g was evaluated from single-crystal EPR data 
using standard formulae 2 ^. The intralayer exchange energy is calculated using Eq. [2] typical uncertainties in the values of J 
resulting from uncertainties in g and B c are ±0.1 K. Neel temperatures In were measured to ±0.04 K using ^iSPe^ 2 -, apart from 
Cu(pyz) 2 (Re04) 2 , where the transition was observed in heat capacity data 2 ^ (typical uncertainty ±0.1 K). The anisotropy 
\J±/J\ calculated using Eq. [3] 



When a sample is within the coil, the signal is V oc 
(dM/di), where t is the time. Numerical integration is 
used to evaluate The sample is mounted within 

a 1.3 mm diameter ampoule that can be moved in and 
out of the coil2£. Accurate values of M are obtained by 
subtracting empty coil data from that measured under 
identical conditions with the sample present. 

Fields were provided by the 65 T short-pulse and 100 T 
multi-shot magnets at NHMFL Los Alamos 2 ^ and a 60 T 
short-pulse magnet at Oxford. The susceptometer was 
placed within a 3 He cryostat providing Ts down to 0.4 K. 
B was measured by integrating the voltage induced in 
a ten-turn coil calibrated by observing the de Haas-van 
Alphcn oscillations of the belly orbits of the copper coils 
of the susceptometer—. 

In the cases that sufficient quantities of material were 
available, Neel temperatures Tn were measured using the 
zero-field muon-spin rotation (/iSR) techniques described 
in Ref32 (see also Refs i 19 ' 24 ); the results are given in 
Table |U Ref.— also discusses the difficulties encountered 
when using more conventional techniques to find Tn in 
quasi-2D magnets. 



III. EXPERIMENTAL RESULTS 

Typical M(H) data are shown in Fig.[l] all compounds 
studied (Table [TJ> behaved in a very similar fashion. At 
higher T, M{B) is convex, showing a gradual approach to 
saturation at high B. However, as T — > 0, the M (B) data 
become concave at lower B, with a sharp, "elbow"-like 
transition to a constant saturation magnetization M sat 
at higher B; no further changes in M occur to fields of 
85 T with the current materials. Wc label the field at 



which the "elbow" occurs B c . As shown in Fig. QJ>, B c 
depends on the crystal's orientation in the field. However, 
in such cases, the M data become identical to within 
experimental accuracy when plotted as Af/M sat versus 
gB, where g is the g-factor appropriate for that direction 
of B (Fig. [He) and Table HJ). This suggests that the g- 
factor anisotropy is responsible for the observed angle 
dependence of B c . 

IV. MONTE-CARLO SIMULATIONS OF 
MAGNETIZATION DATA 

The magnetic properties of the materials in Table U 
are well described by S = \ Cu 2+ spins on a square lat- 
tice. The layers are arranged in a tetragonal structure^; 
coupling between the layers depends on the packing of 
the molecules between. To model data such as those in 
Fig. [TJ we assume that the interaction between the spins 
is purely Heiscnberg-likoi', resulting in the Hamiltonian 

h = j J2 s, -s, • ./. X s -- s j - h J2 s * w 

where J (J±) is the strength of the intra- (inter-) pla- 
nar coupling and h = gi-J-sB is the Zeeman energy pro- 
vided by the (uniform) magnetic induction B. The first 
(second) summation refers to summing over all nearest 
neighbors parallel (perpendicular) to the 2D xy-plane. 
The dependence of M as a function of B is studied for 
a range of values of J± / J [from J± = (completely de- 
coupled layers) to J±/J = 1] using large-scale numerical 
simulations. Note that the orientation of B only affects 
h through the anisotropy of the g-factor, in agreement 
with the discussion of FigQJc) above. 
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FIG. 1: (a) Magnetization M of [Cu(HF 2 )(pyz)2](BF 4 ) pow- 
der versus field B; data for T = 0.5, 1.5 and 4.1 K are 
shown (traces for 0.5 and 1.5 K overlie), (b) Magnetization 
of [CuF2(pyz)](H20)2 single crystals with B applied parallel 
(upper 4 traces) and perpendicular (lower trace, 0.5 K) to the 
2D planes. Data for T = 0.5, 1.5, 4.1 and 10 K are shown 
for the Bii case. The "elbow" denoting saturation occurs at 
Be = 28.8 T for B\\ and B c = 33.1 T for B±. (c) Normal- 
ized M data (T = 0.5 K) for [Cu(HF 2 )(pyz) 2 ]C104 single 
crystals versus gB, where g is the appropriate g-factor, for B 
parallel(B||) and perpendicular (B±) to the layers. 



The stochastic series expansion (SSE) metho d 30 ' 31 is a 
finite-T Quantum Monte Carlo (QMC) technique based 
on importance sampling of the diagonal matrix elements 
of the density matrix e~@ H . Using the "operator- loop" 
cluster update 3 -^, the autocorrelation time for system 
sizes considered here (up to « 3 x 10 4 spins) is at most 
a few Monte Carlo sweeps even at the critical T for the 
onset of magnetic order—. Estimates of ground state ob- 
scrvables are obtained by using sufficiently large values of 
the inverse T, f3. We have further found that the statis- 
tics of the data obtained can be significantly improved by 
the use of a tempering schem o 33 ' 34 . We use parallel tem- 
pering 3 ^, where simulations are run simultaneously on a 
parallel computer, using a fixed value of J_l and different, 
but closely spaced, values of h over the entire range of 
fields up to saturation. Along with the usual Monte Carlo 



updates, we attempt to swap the values of fields for SSE 
configurations (processes) with adjacent values of h at 
regular intervals (typically after every Monte Carlo step, 
each time attempting several hundred swaps) according 
to a scheme that maintains detailed balance in the space 
of the parallel simulations. This has favorable effects on 
the simulation dynamics, and reduces the overall statisti- 
cal errors (at the cost of introducing correlations between 
the errors, of minor significance here). Implementation 
of tempering schemes in the context of the SSE method 
is discussed in Ref<^. 
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B/B c 

FIG. 2: (a) Calculated magnetization M of a spin-^ square 
lattice with added interlayer anisotropy (see Eq. [TJ. Simula- 
tions are shown for J±/J — (uppermost curve) ^, \, \ 
and 1 (lowest curve). Increasing J±/J raises the critical field 
for saturation B c and reduces the curvature of M(B) below 
B c . (b) Comparison of the model shown in (a) (curves) with 
experimental data (points) for a [Cu(HF)2(pyz)2]C104 crys- 
tal (T = 0.5 K, field applied parallel to the 2D planes). Both 
model and data are plotted in reduced units M/M sa t, B/B c . 
Note that the experimental data lie between the J±/J — 
and J±/J= jg model curves. 



V. COMPARISONS OF MODEL AND DATA 

Fig. [DJa) shows the predictions of the model for low 
T, and Fig.[2fb) shows a comparison with typical experi- 
mental data. This is made by plotting both model results 
and experimental data in dimensionlcss units, M/M sa ± 
and B/B c . As M sa t is known, there is in effect only one 
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fit parameter, B c . The value of B c is varied until there 
is a satisfactory overlap of the data and one of the model 
curves^. The curvature of the data and the presence of 
the "elbow" place tight constraints on such a compari- 
son; that in Fig[2jb) is typical of all of the materials in 
Table[H with their M(B) data falling between, or closest 
to, the J±/J = or J±/J = jg numerical curves, indi- 
cating a high degree of anisotropy. We shall give further 
justification for this assertion below. 




i 

50 100 150 200 250 300 
Temperature (K) 




2 4 6 8 10 12 



Exchange energy magnitude \J | (K) 

FIG. 3: (a) Temperature dependence of the low-field sus- 
ceptibility x of [Cu(HF)2(pyz)2]SbF6 (points) fitted to the 
2D Heisenberg expression (curve) of Ref.— to obtain an esti- 
mate of J. (b) Experimental saturation fields times p-factor 
(gB c ) plotted against values of J deduced using fits such as 
that in (a) (points). The straight-line fit to the data yields 
gB c /\J\ = 6.2 ±0.2 TK" 1 . 

For such highly anisotropic magnets, the model pre- 
dicts the ratio gB c /\J\ to take values in the range 
5.95 TFT 1 (J ± = 0) to 6.10 Thr 1 (J ± /J = i). Since 
the experimental uncertainties involved in the location 
of B c are ~ 1 — 2%, and the errors in g are ~ 1%, no 
significant loss in accuracy occurs if we employ the mean 
value, 

^ « 6.03 TFT 1 , (2) 

in what follows. To check the model prediction, a fit 
of the T-depcndcnt low-field susceptibility following the 
method of Refi^i was used to determine J independently 
for a selection of compounds (Fig. [3{a))2£; the values 



obtained are compared with gB c in Fig [^b). As can 
be seen, the points lie close to the line gB c /\J\ = 6.2 ± 
0.2 TK _1 , in good agreement with the predicted value 
(Eq. [5]). As noted above, it is possible to determine the 
value of B c to a very good accuracy (±1 — 2%); once the 
dimensionality of the magnet in question is seen to fall 
within the limits J±/J ~ — using fits such as those 
in Fig. HJb), Eq. [5] almost certainly presents the most 
accurate method for evaluating J. Intralayer exchange 
energies J derived in this way from measured values of g 
and B c are in Tabled] 

Finally, an estimate of the anisotropy \J±/J\ can be 
made using the results of quantum-Monte-Carlo simula- 
tions of quasi-2D Heisenberg antiferromagnets 3 ^, which 
found 

™= ^(2.43-2.30* B), < 3 > 

where both |J| and Tn are measured in K; the result- 
ing values are given in Table fl] Note that Eq. [3] is a 
rapidly- varying function of |</|/Tn; small shifts in cither 
parameter result in quite large changes in | J±/J\. Given 
the experimental and other errors in J and Tn (~ a few 
%), the derived values of \J±/J\ will probably bc within 
a factor ~ 2 of the true values. In spite of this caveat, the 
\J±/J\ values in Tableware ah ~ 0.01, in good agreement 
with the fits to the magnetization data (e.g. Fig [2j) that 
suggested < |Jj_/J| < j§ for all of the compounds. 

VI. SYSTEMATIC TRENDS IN THE 
[CU(HF 2 )(PYZ) 2 ]X FAMILY 

Having established a reliable method for deriving J 

and the anisotropy from M(B) and Tn, we focus our 

remaining discussion on the compounds with formula 
[Cu(HF 2 )(pyz) 2 ]Xi&23 where x can bc CIQ ^ B p 4; pFe 

etc. (Table UJ. All of these materials possess very 
similar extended polymeric structures consisting of 2D, 
four-fold symmetric [Cu(pyz)2] 2+ sheets in the afr-planes 
that are connected along the c-axis by linearly-bridging 
[HF2] - ions (Fig.[H) 18 i 23 . The Cu-Cu separations are sim- 
ilar along the Cu-(pyz)-Cu (0.6852 nm) and Cu-FHF-Cu 
(0.6619 nm) linkages, so that the structure may be de- 
scribed as pseudo-cubi o 18 i 23 ; the X anions are placed in 
the body-centre positions within each "cube" . The Cu-F 
and Cu-(pyz) bonds result in the Cu 2+ d x 2_ y 2 orbitals 
lying within the ab planes, as evidenced by the (/-factor 
anisotropy observed in EPR measurements 2 ^; as noted 
above, the ab planes also correspond to the 2D planes 
within which the strong exchange pathways occur. 

There is little variation of the lattice parameters across 
the [Cu(HF2)(pyz)2]X serie a 18 ' 23 ; given this similarity, it 
is at first sight surprising that the exchange parameters 
J in Table |T] fall into two distinct groups: the compounds 
with tetrahedral anions (X = BF4 and CIO4) possess 
values J ?s — 7 K and those with octahedral anions (X 
= PF 6 , SbF 6 and AsF 6 ) have J w -13 K. Note also 
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FIG. 4: Experimentally-determine d room-temperature 
crystal structures of [Cu(HF2)(pyz)2]X 2D organic magnets. 
In each figure, the Cu-F-H-F-Cu linkages (parallel to c and 
perpendicular to the ac planes) are vertical. The upper figure 
is for X = SbFe whilst the lower shows X = BF4. Note that 
for X = BF4 (lower) , the planes of the pyrazine molecules are 
tilted by 31.6° away from being orthogonal to the ab planes. 
However, for X = SbFe (upper), the planes of the pyrazine 
molecules are perpendicular to the ab planes. (The [SbFe]~ 
and [BF4] - anions have been omitted for clarity; they inhabit 
the "body-center" sites of each approximately cubic unit.) 



that the compounds with tetrahedral anions are more 
anisotropic (\J±/J\ ~ 10~ 3 ) than those with octahedral 
anions (\J±/J\ ~ HT 2 ). 

We first discuss whether the X anions can play a direct 
role as exchange pathways between the Cu 2+ ions. First 
of all, the anions [CIO4] - and [BF 4 ] - have radically dif- 
ferent electronic orbitals, and yet the intralayer exchange 
energies for the magnets containing these anions are very 
similar (see the first three rows of Table HJ. Moreover, 
as mentioned above, the X anions are not within the 
2D Cu 2+ planes, but at body-center positions within the 
"cubes" (Fig. 3]). The Cu-(anion) separation is therefore 
roughly the same for the interlaycr and intralayer direc- 



tions; if the anions played a direct role as exchange path- 
ways, one might expect that the [Cu(HF2)(pyz)2]X com- 
pounds would have \J±/J\ values that were somewhat 
larger (i.e. more isotropic) than the values ~ 10~ 3 — 10 -2 
that are actually observed (see Table UJ. Therefore, in- 
stead of playing a direct role in the exchange, it is more 
likely that it is an anion size effect on the crystal struc- 
ture that is affecting the exchange pathways. 

The structural difference that may well be responsi- 
ble for the factor of two change in J is the configura- 
tion of the pyrazine molecules within the Cu-(pyz)-Cu 
linkages. Fig. [3] (upper) shows [Cu(HF 2 )(pyz) 2 ]SbF 6 , 
one of the systems with octahedral anions; the struc- 
tures of X = AsFg, PFg are very similar. At room tem- 
perature, the larger size of the octahedral anions com- 
pared to the tetrahedral examples forces the pyrazine 
ligands to stand up virtually perpendicular to the ab 
planes^ 3 - (Fig. By contrast, the planes of the pyrazine 
ligands in the compounds with the smaller tetrahe- 
dral anions, [Cu(HF) 2 (pyz) 2 ]BF 4 (Fig. (lower)) and 
[Cu(HF) 2 (pyz) 2 ]ClC>4, are not perpendicular to the ab 
plane, but are tilted away by 31.6° (X=BF 4 ) or 25.8° 
(X=C1C>4) in a pattern that preserves the four-fold sym- 
metry of the Cu 2+ sites^ 3 -. As one lowers the tempera- 
ture 2 ^, the tilt of the pyrazine ligands increases in all of 
the compounds (i.e. those with octahedral and tetrahe- 
dral anions); however, in the case of the octahedral an- 
ions, space-filling considerations restrict the tilt to values 
in the range 5° (X = PF 6 ) to 12° (X = AsF 6 ). This sug- 
gests that it may be the orientation dependence of the 
pyrazine ligand that produces the factor two difference 
in J, with the more perpendicularly-disposed pyrazincs 
(Fig. [5]) presenting a more efficient exchange pathway 
within the layers. 



VII. SUMMARY 

Magnetization experiments have been carried out on 
nine organic Cu-based 2D Heisenberg magnets. These 
systems exhibit a low-T magnetization that is concave as 
a function of field, with a sharp "elbow" transition to a 
constant saturation value at a critical field B c . Monte- 
Carlo simulations including interlayer exchange quanti- 
tatively reproduce the data; the concavity indicates the 
effective dimensionality and B c is an accurate measure 
of the in-plane exchange energy J. Taken in conjunction 
with Ncel temperatures derived from muon-spin rotation, 
the values of J may be used to obtain quantitative esti- 
mates of the exchange anisotropy, | J±/ J\. 

We suggest that in organic magnets of the form 
[Cu(HF 2 )(pyz) 2 ]X, where X is an anion molecule, the 
sizes of J and Jj_ may be controlled by the tilting of the 
pyrazine molecule with respect to the 2D planes. Thus, 
it may be possible to use molecular architecture to design 
magnets with very specific values of J, tailored to a par- 
ticular desired property. One possible application would 
be an organic magnet designed to simulate the antifer- 
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romagnctic interactions germane to cuprate supercon- 
ductivity 3 - ' 11 ^ 11 ' 12 ! 13 ! 14 , but with exchange energies small 
enough to permit manipulation of the magnetic ground- 
state using standard laboratory fields. 
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